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In this paper we link the physics of topological non-linear a models with that of Chern-Simons 
insulators. We show that corresponding to every 2n-dimensional Chern-Simons insulator there is a 
(n — l)-dimensional topological non-linear a model with the Wess-Zumino-Witten term. Breaking 
internal symmetry in these non-linear a models leads to non-linear a models with the 8 term. [This 
is analogous to the dimension reduction leading from 2n-dimensional Chern-Simons insulators to 
(2n — 1) and (2n — 2)-dimensional topological insulators protected by discrete symmetries.] The 
correspondence described in this paper allows one to derive the topological term in a theory in- 
volving fermions and order parameters (we shall referred to them as "fermion-cr models" ) when the 
conventional gradient expansion method fails. We also discuss the quantum number of solitons in 
topological non-linear a model and the electromagnetic action of the (27i — l)-dimensional topological 
insulators. Throughout the paper we use a simple model to illustrate how things work. 



I. INTRODUCTION 

In the past couple of years the topic of topological band 
insulators (TBI) protected by time reversal symmetry^"^ 
has attracted considerable interests. These insulators 
are described by free fcrmion Haniiltonians whose band 
structure are topologically non-trivial. (Because the 
mean-field quasiparticle Hamiltonian of a superconduc- 
tor is also a free fermion theory, the notion of topo- 
logical band insulator can be simply generalized to in- 
clude "topological superconductors"^"^.) Two and three 
dimensional TBI with Dirac-like edge/surface band(s) 
have been prcdicted^""^^ and discovered^^"^^. In addi- 
tion, a classification of TBI in different spatial dimen- 
sions has been achieved^'^'^^. A dimensional reduction 
scheme that allows one to generate lower dimensional 
TBI from even-dimensional Chern-Simons insulators has 
been cstablished^^'^^. 

On a different front, the effects of topological terms 
on the dynamics of Goldstone modes and the quantum 
number of solitons and instantons in non-linear a (NLcr) 
models have a long history, and continue to attract strong 
interests from the physics community^*'. One of the earli- 
est examples is the charge associated with the solitons of 
the Pcicrls-Frohlich order parameter in one-dimensional 
charge (bond) density wave systems^ ^"^'*. This arises 
when the fermion field, which couples to a topologically 
non-trivial order parameter profile, is integrated out. In 
terms of the fermion spectrum this amounts to the pres- 
ence of in-gap eigenstates (or in short zero modes) local- 
ized on the soliton. Another well-known example is the 
0-term in the non-linear a model describing one dimen- 
sional antiferromagnetic spin chains^^. There, depending 
on whether = tt or [Mod(27r)], the half-integer or inte- 
ger spin chains are gapless or gapped. In condensed mat- 
ter physics the appreciation of the physical meaning of 
various topological terms grows as the number of realiza- 
tion of topological NLcr models increases. For example, 
recently it is appreciated that, in a NLcr model describing 
competing orders, one of the effects of the Wess-Zumino- 



Witten (WZW) term is to make the topological defects 
of an order parameter to carry the quantum number of 
its competition order parameter^^"^^. 

In many-body physics one often encounters La- 
grangians which arc quadratic in the fermion fields, and 
in which some fermion bilincars couple to space-time de- 
pendent order parameters. This type of action can be 
viewed as the "fiuctuation-corrected mean-field theory" 
of an interacting fermion problem. In a seminal paper 
Abanov and Wiegmann'^*^ gave a systematic derivation of 
several kinds of topological terms and soliton quantum 
numbers when the fermions in such fermion-cr models 
arc integrated out. Their method is designed for models 
where, without the order parameter, the fermions have a 
gapless Dirac spectrum. The order parameters introduce 
various types of "mass term" (i.e., fermion bi-linears) 
each of which can gap out the Dirac node(s). Under the 
condition where the number of such mass terms is right, 
and when the mass matrices anti-commutc with the Dirac 
matrices and with one another, a topological term (which 
is a pure imaginary term in the Euclidean action) is gen- 
erated upon integrating out the fermions. The topologi- 
cal term can fundamentally affect the ground state prop- 
erties and the low energy dynamics of the order parame- 
ters. However, in condensed matter physics gapless sys- 
tems are not always described by the massless Dirac the- 
ory. For these systems the method presented in Ref. 30 
may not work and a more general approach is desired. 
Moreover, on a conceptual level, it will be very satisfy- 
ing to understand the relation between two of the areas 
in topological condensed matter physics: the topological 
band insulators and the topological NLcr models. This 
paper serves to fulfill these goals. The main result of this 
paper can be summarized as follows. 

Let Hq be a translationally-invariant free fermion 
Hamiltonian defined on a (n— l)-dimensional lattice sys- 
tem. Now let us add fermion bilinears that couple to 
a set of {n + 2) space-time dependent order parameters 
ma{x, t) (a = 1, n + 2), and hence produce the follow- 
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ing Euclidean action: 



S = 



dt 



i 



a=l Lj 



(1) 



where i, j label the lattice sites and x is the averaged po- 
sition of site i and j; ipj = ("0111 V'Id) f-re D fermion 
creation operators on site i. We will focus on local cou- 
pling so that the D x D matrices M,"- vanish sufhciently 
fast as sites i and j get further apart. The question we are 
interested in answering is under what condition a WZW 
term for the ma 's will be generated upon integrating out 
the fermions. In the case where Hq possesses a gapless 
Dirac spectrum the answer is provided by Ref. [30] . More 
explicitly, if the low energy and long wavelength physics 
in the continuum limit is described by the following ac- 
tion 

/n— 1 
dtd'^x ij^'dttp + Y.'PH-^da^"')^ 



n+2 



(2) 



where r"'s arc anti-commuting matrices namely 
{r",r'^} = 25"^ and Af" are constant mass matrices. 
The condition for the generation of the WZW term is 
when the matrices r"'s and the mass matrices A/°'s sat- 
isfy 



_iai---a„+2 



(3) 



Here e is a totally anti-symmetric tensor. 

Gapless but non-Dirac dispersions often appear in con- 
densed matter systems. For instance, Hq can describe 
quadratic or high order band touching, or even possesses 
a Fermi surface. For these cases under what condition 
will a WZW term be generated upon integrating out 
fermions? 

In this paper we shall show that if integrating out the 
fermion field in Eq. (1) generates a level k WZW term 



27rfc 



^wzw 



Area(5'"+i) Jq 
X na^duTia^donasdinai ■ • •5„_ina„+2, (4) 



where Ua = ma/\/X]6 '^6' then there exists a finite range 
of A where the following Bloch Hamiltonian 

H{ki, ■ • • , fc„_i;pi, • ■ • ,Pn+i) = Hoiki, ■■■ , fc„_i) 

n+2 

+X'^ma{pi, - ■ ■ ,_p„+i)M"(fci, ■ • • ,A:„_i) (5) 

a=l 



describe a 2n-dimensional Chern-Simons insulator, so 
long as the Pontryagin indcx'^^ P of the map 



is non-zero. 



In Eq. (5) M^'iki,-' 



,"i„+2(p)) (6) 
,kn-i) is the Fourier 



transform of A/^" . In Eq. (4) an auxiliary coordinate u is 
introduced so that 

n{t,xi, ■ ■ ■ ,x„_i,u = 0) = (0, 0, • • • , 1) 

h{t,xi, ■ ■ ■ = 1) = n{t,xi, ■ ■ ■ ,Xn-i)- (7) 

The Pontryagin index of the map in Eq. (6) is given by 



P = 



^ „ n+l 



Area(S'"+i) 



X Tld^ ^pi ^P2 ^0-3 • ■•^Pn + l 1 



(8) 



and the n-th Chern number characterizing the 2n- 
dimcnsional Chern-Simons insulator is given by-'^^ 



n!2"(27r)' 



d^''k(} 



(9) 



in which fij is the momentum space Berry curvature as- 
sociated with the occupied bands. We emphasize that 
the above result is valid regardless of whether the gap- 
less Hamiltonian Hq is Dirac-like or not. Of course once 
we derived the WZW term, other topological NLtr model 
can be obtained by breaking the internal 0{n + 2) sym- 
metry. In this way a link between the topological NLct 
models and Chern-Simons insulators is established. 

Conversely if Eq. (5) describes a Chern-Simons insu- 
lator for a range of A, it is possible to pick an appropri- 
ately bounded order parameter field rn{t,x), so that the 
fermion NLcr model given in Eq. (1) has a WZW term 
upon fermion integration. The readers are referred to 
appendix A for a proof of the above correspondence. 

Actually, the connection between 2n-dimensional 
Chern-Simons insulators and (n — l)-dimensional topo- 
logical NLct models discussed above is a special case 
of the following more general relation. Consider a 
277,-dimensional Bloch Hamiltonian describing a Chern- 
Simons insulators with a nonzero C„ for a finite range of 
A (0 < A < A*) 
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l) = Hoih 



,kn-l) + V(ki 



; Ami(pi, ...,pn+i), Am„+2(pi, (10) 



If the map (pi, • • • ,Pn+i) (jTii, • • • , m„+2) has a 
nonzero pontryagin index P, then the following steps gen- 
erate a (n— l)-diniensional fermion-cr model which yields 
a level k = C„/P WZW term upon integrating out the 
fermions. (i) Perform the Fourier transform with respect 
to fci, • • • , kn-i to rewrite Eq. (10) in real space, (ii) Re- 
place TTiaipi-,--- ,Pn+i) with an appropriately bounded 
space-time dependent fields ma{t,x). 

The correspondence discussed above establishes a link 
between topological NLct models and Chern-Simons in- 
sulators. It also allows one to determine whether inte- 
grating out fermion field generates topological terms in 
the NLcr model describing the order parameter dynamics 
when the conventional gradient expansion of in Ref. [30] 
fails. The rigorous proof for the above correspondence is 
presented in appendix A. 

In the rest of the paper wc shall use a simple model 
to illustrate the principal ideas. It turns out that for 
this simple model the Abanov-Wiegmann method does 
not work, consequently the scheme we discussed above 
is necessary to derive its WZW term. 



II. A SIMPLE MODEL WITH QUADRATIC 
BAND TOUCHING (QBT) 

As a simple example of a gaplcss Hamiltonian which 
has a dispersion different from that of massless Dirac 
fermions we consider the following one dimensional lat- 
tice model (this model is motivated by the spectrum of 
bilayer graphene). The schematic representation of the 
model is shown in Fig. 1(a). Here back and red indicate 
strong and weak bonds respectively. The vertical blue 
bonds have the strongest hopping. There are four sites 
per unit cell and the 4x4 Bloch Hamiltonian is given by 



by 



H{k) 



/ A Tisik) 

A T2A{k) 

T^sik) 

72*4 (fc) 



(11) 



where ri3(fc) = (1 - St) + (1 + 6t)e-''' and T24,{k) ^ 
{l + 6t) + {l- (5t)e"*'=. Here l±St and A arc the strength 
of black, red and the vertical bonds respectively. When 
dt = and A = 2 the dispersion relation is shown in 
Fig. 1(b). Note that two bands touch quadraticalli/^^'^^ , 
at fc = TT. When A ^ 1 the Hamiltonian in Eq. (11) 
can be reduced to an effective 2x2 Hamiltonian by per- 
forming second order perturbation theory to remove the 
strongly hybridizing vertical bonds. The result is given 




cos 4 - iSt s'm ^] 



(cos I + iSt sin |) 




(12) 



The corresponding dispersion is shown in part (c) of 
Fig. 1. In the following we shall multiply the Hamiltonian 
in Eq. (12) by — ^ and consider the effective model 



Hik) 







(cos I + iSt sin |^ 



\ (cos I — iSt sin ^) 
(1 - 6t^) + il + Sf)cosk 



J 

[25t sink] T2. 

(13) 



One might argue that the QBT point can be viewed 
as two Dirac points merged together. Hence upon inte- 
grating out fermions one can still apply the method of 
Ref. [30] , and sum the results associated with each Dirac 
point at the end. This is in fact not generally true even 
though there are cases the conventional gradient expan- 
sion work^^. As the Dirac point merges, the momentum 
range with linear dispersion vanishes. After the Dirac 
points merge the length scale above which one can think 
of the system as two separated Dirac points diverge. It 
is therefore by no means clear that one can simply sum 
the results assuming the Dirac points are separated to 
determine those for the QBT model. 



III. THE SOLITON OF THE QBT MODEL 

A. The Jackiw-Rebbi soliton 

Setting St ^ in Eq. (13) opens a gap in the dis- 
persion (Fig. 2(a)). In the classic work of Jackiw and 
Rebbi^^ it was shown that in a soliton mass background, 
the one-dimensional Dirac Hamiltonian exhibits one mid- 
gap state per soliton. When we let the gap parameter St 
in Eq. (13) to assume a soliton profile (see Fig. 2(b)) 
mid-gap states appear in the energy spectrum (Fig. 2(c). 
However unlike the Dirac Hamiltonian, there are two mid 
gap states per soliton, as expected for the QBT. It is 
worth to point out that on this specific lattice the mid- 
gap states are true zero modes; they are protected by the 
"chiral" symmetry T^H{k)T^ = —H{k) in Eq. (13). 

B. The zero mode theorem 

In this section we prove the zero mode theorem. In 
order to allow a comparison with the Jackiw-Rebbi zero 




(b) E (c) E 




FIG. 1: (color on-line)(a) A schematic representation of the 
QBT model. Here back and red indicate strong and weak 
bonds respectively. The vertical blue bonds are much stronger 
any other hopping, (b) The dispersion relation for the special 
case when the back and red bonds have equal strength, (c) 
The dispersion of the model in Eq. (13) when the strongly 
hybridizing dimers are eliminated. 
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FIG. 2: (color on-line) (a) In the presence of a non-zero St 
an energy gap opens at the quadratic touching point. Here 
St — 0.2 is used, (b) A soliton 5t{x) profile, (c) Schematic 
representation of the energy spectrum associated with a soli- 
ton St{x) profiole. The red line marks the position of mid-gap 
soliton states, (d) The eigenfunctions associated with the zero 
modes. 



mode thcorcm^^, we work in the continuum limit. By 
expanding Eq. (13) around fc = tt we obtain 



(14) 



where q — k — -k. When 5t is position-dependent the 
above Hamiltonian becomes 



with 







Q = i2St{x) + idx 



(15) 



(16) 



Since the third Pauli matrix, T3, anticommutes with H in 
Eq. (15), the zero energy eigenfunctions must be eigen- 
vectors of it. Thus they are either of the form | )^ ^ ] or 

u{x) 


Let us first consider the first possibility 




(17) 



In this case Q^v{x) = and there are two linearly inde- 
pendent solutions, namely. 



v{x) ^e-^fo^*^y^''y, and 
v{x)^xe-^^o'''-y'>''y. 



(18) 



Note that for the St{x) profile shown in Fig. 2(b) (where 
the center of the soliton is at a; = 0) these are decaying 
solutions as a; ^ ±00, hence are normalizable. Note 
that the first line of Eq. (18) is the same as the Jackiw- 
Rebbi soliton solution. For St{x) = Sto ta.nh{x / £,) these 
two solutions are illustrated as the blue and red curves 
in Fig. 2(d). 

The second possibility 











u{x) 




= 0, 



(19) 



yields the solutions 

Uix) -fo ^t{y)dy ^ 

u{x)^xe^^o^'^y)dy. 



(20) 



For the St{x) profile shown in Fig. 2(b) these two solution 
are not normalizable, hence must be rejected. 



C. The spectral flow associated with the 
Goldstone-Wilczek soliton 

Goldstone and Wilczck'^'^ generalized the work of 
Jackiw and Rebbi^^ and considered a model with two 
spatial-dependent mass terms. Viewing the two mass 
parameters as the components of a 2-vector, they consid- 
ered the situation where the mass vector rotates as the 
position X is varied. In the case where the ID system is 
a closed ring there are soliton configurations where the 
mass vector wind integer number of times as x traverses 
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the ring. In Ref. 30 it is shown that for winding number 
n the sohton charge is \AQ\ = ne. 

In exactly the same fashion we generahze the QBT 
model of Eq. (13) to include one more mass term. In the 
situation where both mass terms (mi, 7712) are constant, 
the Hamiltonian read 



Hik) 



TO2 (cos 2 + sin 2) 1(21) 

- irrii sin ^) — / 



7712 represents the difference of on-site potential between 
the two chains. First, we let (7711,7712) 777(005 0, sin 6*), 
and plot the energy spectrum of H{k, 9) as a function 
of k and 9. The results is shown in Fig. 3(a) where 777 
is chosen to be 0.2. Clearly, for all values of 9 there is 
an energy gap. Next, we consider the case where 9 is 
position-dependent and winds by 27r7i as x traverses the 
ID ring. In Fig. 3(b) We show the spectral flow across 
= as 77 is varied. (For non-integer 77 there are spatial 
discontinuity in the mass terms.) Clearly in every 27r 
cycle of 9, there are two energy levels cross E = from 
above to below. This is in contrast with the Goldstone- 
Wilczek model^'^ where only one level crosses E ^ per 
2tt cycle. 

The above spectral flow is summarized by the following 
term in the imaginary-time effective action 



27r 



X 2 / dxAodJ. 



(22) 



The same term with half of the coefficient, i.e., 
— ^ J dxAodx9 was derived first by Goldstone- 
Wilczek^^, and summarized in Abanov-Wiegmann"^" 
as one of the applications of the gradient expansion 
method. 

Thus, to derive the effective action in Eq. (22), one 
naturally would attempt to generalize the gradient ex- 
pansion method to the continuum version of the model 
studied in this section, namely. 



H 



27772 (x) 
{idx — 727711 (x))^ 



(idx + i2mi{x))'^ 
-2m2{x) 



(23) 



(b) 



0.15 




-0.15 




FIG. 3: (color on-line) (a) The energy spectrum of Eq. (21), 
with mi — >■ 777 cos 7772 — >■ msmO, as a function of k and 6. 
Notice that for all values of 9 there is an energy gap in the 
ID dispersion E{k). (b) The spectral flow as as function of 
the winding number 71. 



couples with the momentum. Nonetheless, we cannot 
rule out the possibility that some variant of the gradient 
expansion method can work for the model in Eq. (23). 
Instead of finding a variant of the gradient expansion 
method, we propose a dimension reduction scheme to 
derive topological terms in fermionic sigma models. The 
dimension reduction scheme is qualitatively different 
from the previous gradient expansion method and is a 
systematic method applicable to fermionic sigma models 
whose i?o(fc) can have generic fermion dispersions such 
as massless Dirac, quadratic touching, or even Fermi 
surfaces. 



IV. A 2D CHERN-SIMONS INSULATOR AND 
THE CORRESPONDING OD TOPOLOGICAL 
NLa MODELS 



D. Beyond the gradient expansion 

It is simple to see that a direct application of the 
gradient expansion method employed previously (e.g. 
in Ref. 30) does not work for the model in Eq. (23). 
The problem is rooted in the facts that (1) the two 
masses 7771 and 7772 do not enter the Hamiltonian on 
equal footing. Because of this, it is not possible to find 
a chiral rotation'^^^'^^ so that locally the mass vector 
is rotated to a reference value. (2) When we expand 
{idx + i2mi{x)) the cross term mixes the spatial 
gradient and rriiix). In other words, the mass terms 



In Ref. 18, Qi, Hughes, and Zhang put forward a 
dimension reduction scheme which allows one to start 
from Chern-Simons insulators in even dimensions and 
obtain lower dimensional topological insulators by 
treating the extra momentum variables as parameters. 
Here we take their program one step further and show 
how to obtain topological fermion-NLa models from 
dimension reduction. We will point out some qualitative 
differences between our view and that of Qi et al as we 
go along. Let us start with the relation between a 2D 
Chern-Simons insulator and the zero-dimensional NLcr 
models as a warm up. 
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A. From the QBT model to a two-dimensional 
Chern-Simons insulator 



kis[0,0.95-],k2MO,2/-7] k, e [0,1.95'-?], k26[0,2/r] 

FIG. 4; (color on-line) Left panel: the image of ki G 
[0, 0.957r], ^2 £ [0, 27r). Right panel: the image of ki G 
[tt, 1.957r], ^2 G [0,27r). We deliberately removed two small 
intervals in ki, namely, (0.957r,7r) for the left and (I.QStt, 27r) 
for the right panel, to reveal the origin (red dot) where H is 
degenerate. 



From Fig. 3(a) we see that by regarding 6 as the 
momentum in an extra dimension, the bandstrcutre is 
that of a 2D insulator. What is less obvious is that 
this insulator is a Chern-Simons (or Hall) insulator with 
axy = or equivalently with the first Chcrn number 
(TKNN index^"^) Ci = 2. This can be shown by directly 
computing the first Chern number 



, f .2, ,,- K|9fc.g(k)|^0(i.k|9fc,g(k)|ui,) 



associated with the insulator described by 



Hik) 



sin fco 



[cosf 



tcos k2 sin -^j " 



[cos ^ + im cos fc2 sin ^ 
— m sin(fc2) 



(25) 



In Eq. (24) i,j = 1,2 and u and v labels the occupied 
and unoccupied bands, respectively. More pictorially we 
can expand the 2x2 Hamiltonian in Eq. (25) as the 
linear combination of three Pauli matrices, and view the 
combination coefficients as the coordinates in . In this 
way we have constructed a mapping from the Brillouin 
zone /ci,2 e [0,27r) to R^. In Fig. 4 we illustrate this 
mapping. The left panel corresponds to the image of 
ki e [0, O.OStt] ,k2 S [0, 27r) and right panel is the image of 
ki G [tt, 1.957r], fc2 G [0, 27r). We deliberately remove two 
small intervals in fci, namely, (0.957r,7r) and (1.957r,27r) 
to reveal the origin (red dot) where H is degenerate. 

Clearly the image is a closed surface, and it wraps 
around the origin twice which implies Ci = 2. If we 
couple the electromagnetic gauge field to the fermions 
described by the Hamiltonian in Eq. (25) and integrate 
out the fermions, the familiar Chern-Simons action will 
be obtained 

Scs, = J d^xdief^'^^A^d^Ax. (26) 

This action summarizes the linear response of the Hall 
insulator. 



B. Explaining the spectral flow of the 
Goldstone-Wilczek soliton 



Let us imagine the Hall insulator described by Eq. (25) 
on a two-torus. We insert a magnetic flux L2(f>{xi) 
through the second hole of the torus and restrict Aq, Ai 
to depend on xi and t only. By substituting A2 = (t>{xi) 
and Aolxi^t), Ai{xi,t) into the Chern-Simons action 
Eq. (26) we obtain 



Scs^ ^ -V^i2 / dxidtAo{xi,t)di(j){xi). (27) 



In obtaining Eq. (27) we have taken into account of the 
fact that the time direction is cyclic. Since none of Aq^i 
and <f> depend on X2 , translation in X2 is a symmetry. As 
the result, k2 is a good quantum number, and the system 
decouples into L2 number of one dimensional subsystems 
each characterized by a different fc2. The tight binding 
Hamiltonian of these one dimensional systems are given 
by 
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HiD{k2) = ^ V'l [2TOsin(fc2 + (j^i) T3 + (1 - cos2(fc2 + 00) n V-^ + XI 



v4+i 



+ imcos(fc2 + t/>i) T2 



(28) 



Here ^\ = (cjj, C2j), and Ti_2,3 are the Pauli matrices. In 
the presence of AQ i{xi,t) each of these one dimensional 
system contributes an imaginary (Berry phase) effective 
action Sid which adds up to Eq. (27). Thus 

XS'i_D[fc2 + cj){xi);Aa{xi,t)] 

i C f 

= -^1^2 / dxidtAa{xi,t)di(l>{xi). (29) 

Eq. (29) imphes 

5'l_D[fc2 + (l){xi);Aa{xi,t)] 

,c(fc2) 



27r 



dxidtAo{xi,t)di(j){xi), (30) 



where 



dk2 
"2^ 



C(fc2) = Ci 



(31) 



Now we show that 0(^2) is acutally independent of fc2. 
Let us consider the special ^2 = case of Eq. (30), and 
set (j){xi) = 9{xi) + k2- Under these condition Eq. (30) 
becomes 

SiD[k2 + 9{xi);Aa{xi,t)\ 

= / dxidtAo{xi,t)die{xi), (32) 
27r J 

Compare the above result with the fc2 7^ case of Eq. (30) 
we conclude c(fc2) — c(0) is independent of ^2- Using 
Eq. (31) we conclude 0(^2) = Ci, Vfc2. Thus 



Sid[4>] Aq\ 



iCi 



dxidtAijdi 



(33) 



Eq. (33) is the desired topological term specifying the 
charge of Goldstone-Wilczek soliton, derived without us- 
ing the gradient expansion. A similar consideration was 
given in Ref.18, where a spatial constant but time de- 
pendent k2 was considered. The authors argued that 
the quantum Hall effect in 2D translates into quantized 
charge pumping as the parameter ^2 in the ID model is 
varied. In this way they derived J = Then by 

charge conservation the authors deduced the charge den- 
sity associated with a spatial dependent to be ^ ^ . 

The above method deriving the spectral flow of the 
Goldstone-Wilczek soliton in ID can be straightfor- 
wardly generalized to higher dimensions. To compute 



the fermion number of a topological soliton in a n- 
dimensional fermion a model with n+l masses rUa (i, x) , 
a ~ we compute the n-th Chern number 

Cn of the corresponding 2n-dimensional insulator ob- 
tained by replacing ma{t,x) with ma{pi, ■ ■ ■ ,Pn), where 
Pi are momenta in the extra dimensions. When the map 
(pi, • • • ,p„) ma(pi, • • • ,pn) has Pontryagin index P, 
the fermion number of the fundamental soliton, i.e., soli- 
ton whose real-space Pontryagin index is equal to one, is 
Cn/P- Since the proof of the above statement is simi- 
lar to the one for the WZW term, we will not present it 
explicitly. 



C. Deducing the existence of edge states in the 
QBT model 

In this section we shall deduce informations concerning 
the edge properties of the QBT model. We start with the 
2D Hall insulator (Eq. (25)) on a torus and adiabatically 
switch on a spatially independent magnetic flux L2(l) in, 
say, the second hole of the torus. Thus A2 (j> &nd, 
like in the last section, we shall allow Ao,i to depend on 
xi and t only. Under these conditions the Chern-Simons 
effective action becomes 

Scs, ^ [ dxidte^f'd^Ap. (34) 

27r J 

On the other hand since the X2 translation symmetry 
remains, k2 is a good quantum number. Hence the action 
in Eq. (34) is the sum of the Berry phase of an assembly 
of ID systems, each characterized by a parameter A:2, i.e.. 



iCi4> 



^L2 f dx^dte'^^da^Ap ^Y.^^D[k2 + ^-^A^]. (35) 



Compare the above result with that obtained by setting 
6 to zero we obtain 



{SiD [k2 + 0; Aa] - SiD [k2 ; A^]) 



iCi 



2tt/L2 



dxidte"'^daAf, 



(36) 



The right hand side of Eq. (36) is a total derivative, and 
vanishes in the absence of boundary. Since the imaginary 
time dimension is cyclic, the only chance for it to be non- 
zero is when the system is open in the xi direction. Under 



8 



that condition Eq. (36) becomes 



k2 



2tt/L2 



dt{AoiLi,t)~Ao{0,t)) (37) 



Perform i6/5Ao{Li,t) and i5/6Ao{Q,t) of the right hand 
side we obtain the extra charge locahzed on the right and 
left ends: 



^AQHght(fc2) = Ci 

fe2 

^AQieft(fe)--Ci 



2tt/L2 ' 



2n/L2 



and 



(38) 



Here A(3iight(fc2) and A(5ioft(fe) are the edge state occu- 
pation (or change of polarization'^*) of the adiabatically 
evolved (^2 k2+(j)) ID system minus that of the orig- 
inal {k2) system. In Fig. 5 we show the energy spectrum 
of the 2D Hall insulator with open boundary condition in 
the xi direction. A vertical cut of this spectrum at fixed 
k2 yields the energy spectrum of the ID model associ- 
ated with parameter k2. In particular, the vertical cut 
at fc2 = gives the energy spectrum of the QBT model 
under open boundary condition. Clearly the QBT model 
possesses zero energy edge states. From the degeneracy 
and wavefunctions of the E ~ eigenstates we deduce 
that there is one edge state localized at each end of the 
system. According to Fig. 5 in addition to /c2 = 0, the 
k2 ~ TT system also possesses zero energy edge states. 
This is not so surprising because k2 = k2 = tt corre- 
sponds to reversing the sign of 6t in the QBT model. By 
calculating the eigenfunctions associated with the in-gap 
band we find that the upward {dE{k2)/dk2 > 0) dispers- 
ing branches near both k2 — and tt are associated with 
the right edge, while the downward dispersing ones are 
associated with the left edge. This is also not surpris- 
ing because a Ci = 2 Hall insulator should have two 
co-moving chiral edge branches at each edge. 

In Fig. 6 we zoom in on the in-gap band of Fig. 5. Here 
the blue and red dots represent occupied and empty ^2 
(= 27rn/L2) states respectively. Panel (a) is the occupa- 
tion before the adiabatic evolution and panel (b) is that 
after the adiabatic evolution. For the upward dispersing 
branch k2 ^ k2 + (j> induces an extra number (= 2J/L2 ^ 
of occupied k2 values near each Dirac point (fc2 = 0, tt). 
Similarly, for the downward dispersing branch there are 
exactly the same number of extra empty fc2 values. As 
a result, the Ci in Eq. (38) counts the number of Dirac 
points in Fig. 5. Equivalently it counts the number of 
times the edge states cross Fermi energy {E = 0) as the 
parameter k2 of the ID Hamiltonian sweeps from to 
2tt. (Note that a crossing from below E ^ to above is 
counted as -1-1 while the reverse is counted as —1.) Of 
course, the fact that when a parameter of a ID system 
is varied, charges flow from one end of the system to the 
other is the charge pumping effect^*'^^. 







- 1 








-o.a 


0.5 ^ 


1.0 


^ 1.5 











FIG. 5: (color on-line) The energy spectrum of the 2D Hall 
insulator with open boundaries in the Xi direction. Vertical 
cut of this spectrum at fixed k2 yields the energy spectrum of 
the ID model associated with parameter fc2. The vertical cut 
at ^2 = gives the energy spectrum of the open-boundary 
QBT model. 



(b) k2->-k2 + C 

E 




k2/^ 



FIG. 6: (color on-Hne) (a) The occupation of the edge states. 
Blue dots represent occupied states and the red dots mark 
the unoccupied states, (b) The occupation after the adiabatic 
change ^2 — ^ ^2 + 4>. 



D. The WZW term in d = 



This time we start from the 2D Hall insulator on a 
torus, and insert fluxes Li6i{t) and £2^*2 (^) through both 
holes and fix the gauge so that Aq — 0, Ai 0i{t), and 
A2 02{t). Substitute these changes into the Chern- 
Simons action (Eq. (26)) we obtain 



5, 



CSi 



iCi 



L1L2 / dt e'^'eaiWhit), (39) 



where a,b = 1, 2. Now let us introduce an auxiliary co- 
ordinate u £ [0, 1] so that 9i^2{t,u) smoothly interpolate 
between 9i^2{t,u = 0) = 0, and 9i^2{t,u = 1) = Oi,2{t). 
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0<k,<;r; 0<k2<;r ;r<ki<2/T; 0<k2<'r 



FIG. 7: (color on-line) (a) The map {ki,k2) h{k\,k2) 
where k\^2 £ [0, 27r). 



Then, Eq. (26) can be rewritten as 

Scsi 

= -^LiL2 j dtj^ du e'"'e'''du[ea{t,u)d,eb{t,u)y 



47r 



/ du / dt e^-'e''''df,ea{t,u)dM,u), 



(40) 



where /i = labels the time and n = 1 labels the u co- 
ordinate. In the derivation above, we employed the fact 
that the only boundary is in the u direction so that inte- 



J 



gral of a total derivative can be turned into a boundary 
integral along the w-dircction. 

Since 9i and 62 do not depend on the spatial coor- 
dinates xi and X2, ki and k2 remain as good quantum 
numbers. Thus Eq. (39) should be the sum of the Berry's 
phase generated by L1L2 number of zero-dimensional sys- 
tems. Each of this system is parameterized by fci and ^2 
and is governed by the following time-dependent Hamil- 
tonian 

H{t;ki^2)^B{t,ki,2)-T, (41) 
where the "effective" magnetic field B{t, fci,2) is given by 



Bxit,ki^2) = l-m^cos'^{k2 



+ [1 + cos2(fc2 + 6'2(t))] cos(fci + 01 (i)), 
By{tMa) = -2msin(fci +6'i(t))cos(fc2 + 6'2(t)), 
Br.(tM,2) - 2msin(fc2 +6'2(t)). (42) 

The three components of the effective magnetic field 
are just the three masses of the zero-dimensional sys- 
tems labeled by fci_2- The sum of the Berry's phase 
SB\k\ + ^2 + ^2(^)1 associated with each B(t\ kx^q) 
is equal to ScS\ : namely 

Scs. = E + ^1 ' ^2 + 02 {t)\ ■ (43) 

fcl,fc2 

Because the following identity which computes the first 
Chern number^° Ci as well as the Pontryagin index P of 
the map (^1,^2) n = B{0i,2{t, u)+fci,2)/|S(6'i,2(i, u) + 
ki,2)\ 



Ci=P=^ [ d2fce-6=n„(fci.2 + 0i,2(t, u))^7n,{kx.2 + ei,2(i, m)) J-n,(fci.2 + 0i,2(i, u)), 
47r J oki ok2 



(44) 



is true regardless of the values oi9i{t,u) and 92{t,u) (see Fig. 7), we will show that Ssiki + 9i{t,u),k2 +92{t,u)] is 
a WZW term of a spin-1/2 system. To do so, we rewrite Scsi as follows: 



ScSi 



L1L2 dt du e>"'e"'^df,9ait,u)d,,9f3{t,u), 



iCi 
47r 

-^^j^LiL2 I dt I du I d^k e''''e''f'[d^9^{t,u)d,9p{t,u)] 



e'""fia(fci,2 + 0i.2it, u))dk,nb{ki,2 + 6'i,2(i, u))dk^n^{ki^2 + Oiait, u)) , (45) 



-i^i^LiL2 j dt du I d^k e^''e''f'[d^,9^{t,u)d,9p{t,u)] 



abc 



na{ki^2 + 9i^2{t,u))dg^rn,{ki^2 + 6'i,2(^, w))9e2nc(fci,2 + 6'i,2(i, w)) , (46) 



E 

k\ .k2 



C IP f 

-i-^ I dt j du e''''=na(fci,2 + 6li,2(t, u))dtnb{ki^2 + 6li,2(t, u))dun^{ki^2 + 6'i,2(t, u)) 



(47) 
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where we employed the chain rule of differentiation from 
Eq. (46) to Eq. (47). Because all zero-dimensional sub- 
systems can be adiabatically connected without clos- 
ing the gap, each of them is then characterized by the 
same WZW term with the same level. Since Scsi = 
Efci.fc, SB[ki + 01 (t), fc2 + 02it)] and Ci/P = 1, we con- 
clude that 

SB[n{t)] = -4 /^*/ e'''''nadtmdun,, (48) 

which is nothing but the well-known Berry phase or 
WZW term describing a spin- 1/2 in a time-dependent 
magnetic field, as expected ! 

E. The e term in d = 

The last subsection establishes the fact that the Berry 
phase (WZW term) associated with the Hamiltonian 
H{t) = h{t) ■ T is 

Ss^-'-jdtj^ du e'^^'^nadtUbdun,, (49) 

where h{t, 0) = (1, 0, 0) and n{t, 1) = h{t). Physically Sb 
measures half of the solid angle sustained by h{t,u) at 
u = 1 (i.e. the physical n trajectory) and the north pole 
{n{t, u = 0)). For the case where 

n{t, u = 1) = (sin cos (/)(t), sin sin0(t), cos 6q) (50) 

Eq. (49) becomes 

5..-.(i^)/..|. ,M, 

This is the 6'-term. 

The above is a simple example (n = 1) of the corre- 
spondence between a 2n-dimensional Chern-Simons in- 
sulator and [n — l)-dimensional topological NLcr models. 
This serves as a warm up exercise for the following dis- 
cussions of the correspondence between a 4D (n = 2) 
Chern-Simons insulator and NLct models in ID. 



V. A 4D CHERN-SIMONS INSULATOR AND 
THE CORRESPONDING ID TOPOLOGICAL 
NLcr MODELS 

In a one dimensional chain (whose gapless excitations 
near the right and left Fermi points are described by 
a masslcss Dirac theory) the dimerization and the an- 
tiferromagnetic Neel order parameter form a four vector. 
The quantum NLcr model describing this "super- vector" 
possesses a topological, WZW, term^^. Among others, 
this means the defects of one order carry the quantum 
number of the other order. Indeed, in the presence of 
spin the dimerization domain wall has two (one for each 



spin species) zero modes. Among the four ways of oc- 
cupying them, two lead to charge- neutral spin 1/2 soli- 
tons, and the other two lead to charge ±1 spin zero 
solitons. The condensation of the former (which is fa- 
vored by repulsive interactions) leads to the (algebraic) 
antiferromagnetic order and the condensation of the lat- 
ter (which is favored by attractive interactions) leads to 
(algebraic) charge density wave/superconducting order. 
The former scenario is encapsulated by the dimcrization- 
antifcrromagnetic WZW term discussed above. 

In the presence of spin, the QBT model described 
in section I, has four zero modes per dimerization 
domain wall (two for each spin). There are totally 16 
possible ways to occupy them. Two give rise to charge 
±2 spin zero, eight lead to charge ±1 spin 1/2, three 
give charge spin 1, and another three give charge 
spin solitons. Therefore in this case the destruction 
of the dimerization order by condescension of domain 
walls can lead to a large number of possible states. In 
this section we will focus on the condensation of spin-1 
charge neutral solitons. This should lead to spin-1 
antiferromagnetic state which has a Haldane gap. The 
purpose of the following subsections is to study the 
WZW term associated with the transition between the 
dimerizcd phase and the Haldane phase. In addition, we 
shall also discuss the edge states of the spin-1 chain, the 
6 term and their relation to the physics of 3D TBIs. 



A. The QBT model and a 4D Chern-Simons 
insulator 

In the following we shall find that the d ~1 topological 
NLcr models discussed above is related to a Chern-Simons 
(or quantum Hall) insulator in d = A dimensions, which 
is constructed in Ref. [41]. This is the n = 2 example of 
the correspondence discussed in the introductory section. 
The Hamiltonian for the 4D Chern-Simons insulator is 
defined as follow: 

H{ki,k2,k3,k4) = -[(l-7n^) + (l + m^)cos(fc4)]r5 
— TTiiFi — 7112^2 — TO3F3 — 27714 sin(fc4)F4, (52) 

where 



mi 


= sinfci. 


7712 


= sin ^2 , 


70-3 


= sin fca , 


70,4 


= 777, + cos fci -|- cos ^2 



+ cosfc3, (53) 



and 

ri.2,3 ^Ti(g) cri^2,3, F4 = r2 (8) /, F5 = T3 (g) / (54) 

Here / is 2 x 2 the identity matrix. The Pauli matrices 
T and cr have sublattices (orbitals) and spin degrees of 
freedom respectively. By setting fci^2.3 to zero one can see 
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that the Hamiltonian in Eq. (52) and Eq. (53) reduces 
to the spinful version of Eq. (13), after a 90° rotation 
around the T2 axis. For, e.g., m — —2.5, straightforward 
computation shows that the above Hamihonian has a 
non-zero second Chern number (also see Ref. 18), C2 = 
2, where 



C2 



1 

4^ 



d^kTr [A2/34 - /13/24 + /14/23] (55) 



with the momentum space Berry curvature given by 
^-(k) = ,,..2 (56) 



X ( Pg {k).dk^ H{k).PE (k) .dk^ i/(k) .Pg {k)-fi^iy 



Here we note that because the Hamiltonian is con- 
structed out of r matrices, its eigenvalues have the form 
(-£;, -E, +E, +E). In Eq. (57) Eq is the energy of the 
two occupied levels and Ee is that of two unoccupied 
levels, and the projection operators Pq^Pe are defined 
as follows 



-fG(k) = ^ IWakX^Qkl 
CK— occupied 

Ps(k) = ^ |u^k>(w/3k|- 

/3— empty 



(57) 



After integrating out the fermions, the effective gauge 
action associated with such a Chern-Simons insulator is^® 



'CS2 



, C2 
'247r- 



(58) 



Eq. (52) and Eq. (58) play the roles of Eq. (25) and 
Eq. (26) in the following discussions. 



B. The existence of surface states in 3D TBIs 

We start with the Chern-Simons insulator described 
by Eq. (52) and Eq. (53) on a real-space 4-dimensional 
lattice with periodic boundary conditions (hence form a 
4D torus). Let us insert flux Ln(t) in the fourth hole (hence 
Ai — !• (p) and restrict ^0,1,2,3 to depends on only t and 
2^1,2,3 • Under this condition Eq. (58) becomes 



•S'off = -i 



L4 I d^xdte^^P^d^AxdpA^. 



(59) 



Since the translation symmetry in X4 remains, is a 
good quantum number. Thus the above effective action 
is the sum of the Berry phase action of an assembly of 
3D systems, each labeled by a different /C4: 



S3D [ki + (/>; A^ 



.Co 



87r2 



L4 / d-^xdW^P^d^AxdnA. 



(60) 



Subtract the (f) = 
result, we obtain 



version of Eq. (60) from the above 



{SsD [ki + (f>; A^] - S3D [ki-, A^]) 



= -*§yi4 / d''xdte''^P^d,AxdpA„. (61) 

For non-compact gauge field, since the integrand of the 
right hand side is a pure derivative, it is only nonzero in 
the presence of boundary. With spatial boundaries 



L4J d^xdie^^P^d^AxOpAa 



= -^§7^4 ^ d^xdt n^e^^P" AxdpA,] (62) 

where n^, is the J^th component of the unit surface normal. 
Note that the right hand side of the above expression 
is the first Chern-Simons action defined on the surface 
suggesting surface Hall effect as discussed in Ref. [18]. 
Combining Eq. (61) and Eq. (62) we conclude 



E 

k4 



surface 



(ki 



C2 Li4> 



(63) 



Here /^.Uxy. suifacG(fc4) is the difference in the surface Hall 
conductance between the adiabatically evolved {k^ — s- 
kii + (p) system and the original (^4) system. Eq. (63) 
is in close analogy of Eq. (38) in earlier discussions. 

Since for two dimensional Dirac electrons a change of 
mass from negative to positive induces lS.<Jxy = l/27r, 
the above result strongly suggests that as k^ — > kj^ + (j), 
C2 number of massless surface Dirac Hamiltonians are 
encountered. Let us label the ^4 values at which the 
surface Dirac Hamiltonian is massless by fc^s. As these 
special fc|s are crossed under k^ k^ + (j), the mass of 
the surface Hamiltonians change sign. Eq. (63) requires 
when we compare the masses of the adiabatic evolved 
surface Hamiltonians with the original ones we will find 
C2 X Li4)/2TT mass reversals. 

For the model discussed in Eq. (52) this is indeed true. 
Under the condition 2:3 direction is open, while xi^2 di- 
rections have periodic boundary conditions, the surface 
Hamiltonians read 



^^surf (k) 

V (-fcicr2 



fc20"l) 

V (fcio-2 + ^20-1) ^ 



h ■m{k4)(T3 
m{ki)(T'i 



top surface 
bottom surface 



(64) 



where w is a velocity parameter, fci.2 are the components 
of the surface wavevector. The mass term 7/1(^4) van- 
ishes when ki = k\. Surprisingly, despite the fact that 
the Hamiltonian in Eq. (52) is time-reversal invariant, 
the fc|s are not or tt. Instead, we find a pair of fc| i and 
^4 2 = 27r — fc| ]^ where miki) vanishes. The value of these 
two fc|s value depend on the parameter m in Eq. (53). 
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FIG. 8: (color on-line) The circle represents k4 £ [0, 2n) and 
the red points mark fcl i and fcj ^. The green arcs in the 
vicinity of ^4 = and = n mark the regime in ki where 
the surface bands merge into the bulk continuum. In these 
^4 intervals the surface states on the top and bottom surfaces 
re-hybridize, and the surface Hall conductance vanishes. 



For example for m = -2.5 we find A;|_i = 0.69922827r 
and fc|_2 = 27r - 0.69922827r. In Fig. 9 we show the k^- 
dependence of the energy spectrum of H in Eq. (52) 
under open boundary condition in 3:3, and with (^1,^2) 
fixed at (0,0). In addition to the gapless points fc| ^ and 
fc| 2 there is another important feature worthy of atten- 
tion. In the vicinity of ^4 = and = it (marked by 
green) the surface bands merge into the bulk continuum. 
As a result, the surface states on the top and bottom sur- 
faces re-hybridize, resulting in the vanishing of the sur- 
face Hall conductance. This is illustrated in Fig. 8 where 
the circle represents k^ g [0, 27r) and the red points mark 
fc| I and fc| 2 . The green arcs represent the intervals of k^ 
in which the surface bands merge into the bulk contin- 
uum. It is important to note that there is no gap closing 
at cither ends of these arcs. 

In Fig. 10, we plot the energy spectrum associated with 
kli as a function of ki^2- Note that despite the fact 
that for fc4 ^ or TT time-reversal is broken, the surface 
Dirac point is located at k = (0,0). In addition, at k^ = 
fc| I and fc| 2 the low energy surface Hamiltonians are 
completely time-reversal invariant! The fact there exists 
two fc|s at which k^ ^ ki + (j) induces m{k4) to change 
from negative to positive, is consistent with the fact that 
C2 = 2 in Eq. (63). Because of the fc|s are not time 
reversal invariant, there is no time reversal symmetry 
that protects the gapless Dirac surface states. However, 
finding symmetry protected gapless Dirac surface states 
is not the purpose of this paper. 

C. The WZW action in 1 + 1 D 

The purpose of this section is to derive the NLcr model 
describing the competing dimerization and magnetic or- 



FIG. 9: (color on-line) The spectrum of Eq. (52) at surface 
momenta fci = ^2 = as a function of ^4. The two points at 
which the surface bandgap closes are kl^i and fcl ^. . 



E 




FIG. 10: (color on-line) The surface band structure obtained 
by setting k4 in Eq. (52) to one of the special values (marked 
by the red dots) in Fig. 8. 



der described at the beginning of section IV. 

We start with the 4D topological insulator described 
by Eq. (52) and Eq. (53) and insert magnetic fiuxes 
ii, 2. 3^*1, 2, 3 in the first to the third holes of the 4D torus. 
Moreover, we shall restrict ^1,2,3 to depend only on X4 
and t. By substituting ^1,2,3 ^1,2,3(3^4, i) and fix the 
gauge so that Ao,4 — > 0, Eq. (58) becomes 

Scs, ^ -i-^L^L^L^ J dtdx^e^^'^ead^e^die,, (65) 

where a,b,c — 1,2,3. Since ki^2,3 are good quantum 
numbers, Eq. (65) is the sum of the Berry phase action 
of an assembly of one dimensional systems, each param- 
eterized by {ki,k2,k^), i.e., 

5'id[A:i,2,3 + ^'1,2,3] 

fcl,fe2,fe3 

= ~i-^LiL2L3 I dtdxie-'^'eadoebdiO, 
12tt^ J 

(66) 
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Again, we introduce an auxiliary coordinate u G 
[0,1] so that 6*1, 2. 3(^1 u) smoothly interpolate between 
«'i,2,3(i,w = 0) = and 6'i,2,3(i,w = 1) = 6*1,2,3(0: and 
consider the following action 



C2 



247r2 



TTl, j dtdxi [ due''''''e^"'^dj6adu6bdx6. 



where /i = 0, 1,2 labels the t,X4 and w coordinates, re- 
spectively. Because the integrand of Eq. (67) is a total 
derivative it is equal to a boundary integral. Since bound- 
aries only exist in the u direction we can easily show that 
Eq. (67) reduces to Eq. (65). As the result Eq. (66) im- 
plies 

J1D[Ki,2,3 + fl,2,3j 



(2^) 
C2 



247r2 



dtdxA 



If we use j to label the lattice sites in the X/^ direction, 
the ID Hamiltonian parameterized by fci_2,3 is given by 



H = - 



3 

^V][(i-"^4)r5 + ^m„rc, 



Q = l 



Fs — imiV^, il^j + h.c 



(69) 



where mi^2,3,4 are given by 



TOi = sin(A:i -I- 9i), 
TO2 = sin(A:2 + ^2), 
7713 = sin(A;3 + 6*3), 

TO4 = (m -I- cos(/ci -I- 6*1) -I- cos(/s2 + 6*2) + cos(fc3 + 6I3)) . 

(70) 



In the following we shall set the parameter m to —2.5. 
For ?Tii 2.3,4 = the dispersion of the one dimensional 
Hamiltonian is shown in Fig. 1(c) where both bands are 
doubly degenerate. The question is what kind of Berry 
phase action for m^'s will be generated upon integrat- 
ing out the fermions. (Due to the quadratic touching 
between the conduction and valence bands, the gradient 
expansion method discussed in Ref. [30] no longer works.) 

It can be checked that the following three-torus to 
mapping (^1,^2,^:3) na{ki^2,3, +6*1,2,3) = "1^(^:1,2,3 + 



(67) 6*1,2,3)/^/ I]q=i "^0(^1,2,3 + 6*1,2,3) has Pontryagin index 
P = 1 regardless of the 6*1,2,3 values, where 



P=7^f t\dKe"^'''''nadk,ni,dk,nA,nd. (71) 
-^0 „^i 

Given this fact and Eq. (67), we now show that the de- 
sired Berry phase action generated by integrating out the 
fermion is the WZW action 



Swzw[n{^4,t)] 
1C2 C 



'-§^J^du J dx4dte''''^e-''-<'nad^n,d,n,dxnd. 



(72) 



The proof is similar to what we did earlier in section III, 
namely, just using chain rule of differentiation. If inte- 
grating out fermion indeed generate Eq. (72), consistency 
with Eq. (68) would require 



/ n $ / dx4dte'^''\'''>'''nad^nhd,n,dxnd 

= ~^^ J dtdxij^ due'''''ef'''^d^ead,0bdx0,. (73) 



Since each of the Ua depends on A:i,2,3 -I- 6*1,2,3(^1 ^4, u) we 
apply chain rule to the left hand side of Eq. (73): 



-zC2 f dhdk^dh I dx^dtef^'^^'^'-^nad^ntd^n^dxrid 



Gn J (27r)3 
—iC2 f dkidk2dk^ 



6tt j (27r)3 
-iC2 f dkidk2dk^ 



67: J (2^)- JO 
-iC2 



du j dx4dte'"'^e'"""^nad0^nbd0^ncde^ndid^6ad^epdxey) 
du [ dx4dte'"'^e''''"^nade,nbdg,ncdg,nde'"^'^ {dMi^OpOxO-y) 



du / dx4_dte^ 



-A^a/J^l^ f dk,dk2dkse-''-''nA,nbdk,nA,nd) {8,9^8,9^8x6^) 





^ ' du I dx4dte'"'^e°'^'' {8,9^,8,9 p8x9-y) . 



247r2 







(74) 
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Indeed, it agrees with the right hand side of Eq. (73)! 
The non-linear a model with the C2 = 2 WZW term 

S* = y dx^dt^^lna^^lna + Swzw\na{xi^,t)\ (75) 

describes the critical point between the dimerised phase 
((mo) 7^ 0) and the Haldane phase. This critical the- 
ory is realized in the spin-1 chain with quadratic and 
bi-quadratic exchanges: 

i/ = J^[S,.S,+i-(S,.S,+i)2] . (76) 

i 

It has been exactly solved by Takhtajan^^ and 
Babujan''^'^, and the critical theory is that of level-2 
SU (2) WZW model whose NLcr model form is given by 
Eq. (75). 




50 100 150 200 



FIG. 11: (color on-line) The envelope of the edge wavefunc- 
tions. 



D. The S = \ spin chain, edge states, and the Q 
term in the NLcr model 



In this section we set toi^2.3 ^ cos0o"ii,2.3 and 
sin 00- We consider the following fcrmion a model 



H = -^V][(l-sin'^o)r5 + cos0o^mfcri 



^3 



k=l 



E 



I r 1 + sin^ 60 



Ts — i sin 6'or4 ipj + h.c 



(77) 



where the dimerization order parameter has been set to 
sin 00 and the magnetic order parameters to cos^o™- 
In the absence of the m^'s Eq. (77) describes a 
one-dimensional TBI whose bandstructure is shown in 
Fig. 2(a) (each band is now two-fold degenerate). Under 
the open boundary condition there is a doubly degen- 
erate zero-energy edge states at each end of the chain. 
Using the representation 

ri,2,3 = n cri,2,3, r4 = T2 (8) /, Ts = T3 ® /, (78) 

the edge states can be made eigenstates of ti and CT3. The 
wavefunctions of these edge states are shown in Fig. 11. 
Notice that we have plotted (— l)"")/'™ ('t- is the site index) 
rather than simply ipn- This is because the minimum 
gap occurs at = tt hence the these wavefunctions have 
a fast spatial modulation. In the presence of non-zero 
m = (rhi, 7712-, rh^) the degenerate edge states are split. 
The occupied state has magnetization along the local m 
direction. In the presence of space-time dependent m, 
integrating out the fermions produces space-time depen- 
dent edge magnetic moments each contributing a Berry 
phase of the form given by Eq. (48). This is precisely the 
well-known spin- 1/2 at the ends of the spin-1 chain. 



By letting fi in Eq. (72) smoothly interpolate be- 
tween n{t^X4,u = 1) = (0,0,0,1) and h{t,X4,u = 0) = 
(sm6Q,cos6om{x4,t)) we obtain the 0-action for gapped 
spin chain 



If 

— / dxAdtrh ■ dfiJi x dr.rh 
J " 



(79) 



where 6* = C2[7r — sin(26'o) — 26*0]. When the system has 
some sort of discrete symmetry, such as inversion, which 
requires 1114 = or Oq =0, the coefficient = C^t^- If 
C2 is an odd integer, the corresponding NL-cr model 
is gapless; if C2 is an even integer, the NLcr model is 
gapped. When there is no such discrete symmetry, the 
bare coefficient of the 9 term, in general, deviates from 
odd integral multiple of tt, at low energy it renormalizes 
to the closest fixed point characterized by the term 
coefficient 2n7r, which corresponds to a gapped system. 
In our particular case the relevant value is 27r. This is 
in analogy with the Z2 classification of 3D TBIs and the 
discrete symmetry, such as inversion, which protect the 
gapless modes for = (2n 4- l)7r, plays the role of time 
reversal symmetry in 3D TBIs. 



VI. SUMMARY 

In this paper we have tackled the issue of what is 
the relation between the following two areas of topo- 
logical condensed matter physics: topological band 
insulators and topological non-linear a (NLcr) mod- 
els. We have established a correspondence between a 
2rt-dimensional Chern-Simons band insulators and the 
(n — l)-dimensional NLcr models with the Wcss-Zumino- 
Witten (WZW) term. The correspondence does not de- 
pend on whether in the absence of the order parame- 
ter (the cr fields) the original gapless n — 1 dimensional 
fermion Hamiltonian is Dirac-like or not. We have ap- 
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plied this correspondence to a simple ID model where the 
gapless Hamiltonian exhibits a quadratic band touching. 
In the presence of the dimerization and antiferromagnetic 
order parameters, the fermion integration yields an 0(4) 
non-linear a model whose WZW term has a doubled co- 
efficient. This action describes the critical point between 
the dimerized phase and the Haldane phase. We demon- 
strate that the the edge states of the 5 = 1 spin chain 
are nicely captured if one starts with the edge state of 
the dimerized ID topological band insulator. 
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Appendix A: The proof of the correspondence between a 2n-dimensional Chern-Simons insulator and a 

(n — l)-dimensional topological NLa model 



In this appendix we provide the proof of the Chern-Simons insulators o topological NLct model correspondence 
discussed in the introduction. 

We consider a (n — l)-dimensional fermion-cr model described by the time-dependent Hamiltonian i/(fc; m(i, x)), 
where x = {xi, • • • , Xn-i) and k = (fci, • • • , fcn-i) are the position and momentum in (n — l)-dimensions, and TO(t, x) 
is the space-time dependent [n + 2)-component order-parameter fields with a fixed |TO(t,a;)|. We assume that the 
space and time derivatives of rh(t,x) do not enter the Hamiltonian. When the order-parameter fields rn{t,x) are set 
to zero, the corresponding Hamiltonian Ho{k) = H{k]rh{t,x) — > 0) is gapless. One familiar case is when Ho{k) has 
linear dispersion at low energy. However, our following results apply for more general situations. 

The condition we need to impose to Ho{k) is that for constant TOq the Hamiltonian H[k]m) has a completely 
gapped spectrum, and for < |m| < m* the energy gap does not close. The conditions we impose above apply 
for nearly all physically interesting situations. (Note that, in a continuum model where the linear Dirac dispersion 
extends to -E = ±cxd, to* can be taken to oo since the energy gap never closes for any non-zero |to|. This may not be 
true for lattice fermion models where the bandwidth of i?o(^) is finite.) Physically to* should be of the order of the 
band width of Ho{k). For such H{k; rn{t, x)), we take the following steps to determine whether a WZW term will be 
generated by fermion integration. 

Step 1, we create a 2ri,-dimensional Bloch Hamiltonian H{k;Xm{p))^ where p ~ (pi,-- - ,pn+i) and {k,p) is the 
2n-component momentum, by simply replacing constant to^ in H{k]rn) with Xma{p). Most importantly we choose 
rn{p) such that it has a nonzero Pontryagin index P (Eq. (8)). For instance, the following map 

rrii (p) = sin pi , for i = l,---,n + l, 

n+l 

mn+2{p) = (to - ^ COS Pi) (Al) 

1=1 



has the Pontryagin index (Eq. (8)) equal to one when the parameter to satisfies n~l<m<n+l. Note that 
the parameter A is introduced so that |Am(p)| < to*. For such |Ato(p)|, the Hamiltonian H [k] Xm(j))) describes an 
insulator so that we can go step 2 below. 

Step 2, we check whether the 2n-dimensional Bloch Hamiltonian Hik; Xrn{p)) possesses a non-zero n-th Chern 
number C„ (Eq. (9)). If C„ 0, a level-fc = Cn/P WZW term (Eq. (4)) will be generated upon integrating out 
fermions in the (n — l)-dimensional Hamiltonian H{k]rn{t^x)), so long as |TO(t,i;)| < to*. The following is the proof. 

The electromagnetic gauge action of a 27i-dimensional Chern-Simons insulator contains the following topological 
termi8'44 



(71-1- l)!(27r)" 



dtd^'^x 



n 4 f) A . ■ ■ rl A 



(A2) 



where = 0, 1, • • • ,2n labels t, xi, - ■ ■ , X2„ respectively, e is the {2n + l)d antisymmetric tensor. Now, we choose 



A^(i, xi, ■ • ■ ,X2n) = 0, for ^ = 0, • • • , n - 1; 

Ai+n~lit,Xi, - ■ ■ ,X2n) = 6'j(t,Xi,--- ,Xn-l), for i 



,n + 1. 



(A3) 
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For this choice of gauge potential, the Chern-Simons action becomes 



Scs„ -i- 



n + l)!(27r)'^ 



(A4) 



where = 0, 1, • ■ • , n — 1 represents t, xi, ■ ■ ■ , Xn-i respectively, &i = 1, ■ • • ,n + l, and d"'~^x = dxi ■ ■ ■ dxn-i- (Note 
that a overall n-dependent sign is ignored above and hereafter since it can be absorbed by relabeling the components 
ofn.) We introduce an auxiliary coordinate u G [0, 1] so that 9i(t,x,u) (i = 1, • • • , n + 1) smoothly interpolate between 
Oi{t, x,u = 0) = and 9i{t, a?, u = 1) = 9i{t, x). [x = (xi, • • • , x„_i)]. In the following a; = 0, 1, • • • , n will correspond 
to xi, • ■ • , Xn-i, u respectively, and eOi - n _ ^ Under this condition Eq. (A4) can be rewritten as 



cs„ 



(n + l)!(27r)'- 
+ l)!(27r)'- 

Cn 

+ l)!(27r)'- 



" 2n 
_i—n 


J dtd^-^x j 


" 2n 


j dtd^-\ j 


_i—n 
" 2n 

_i—n 


j dtd"-^x J 



du 



ai ■ ■ -a^iCiri+i J^i---h. 



lg«l-«„+l5^^^^ et,^da,Ob,---dajK^, , (A5) 



ai ■ ■ - aTian-i-i 6i ---fcri 



da„ + iObidai0b2 ■ ■ ■ da„0b„ + i , (A6) 
dai 0bi da2 ^ba ■ ■ ■ f^a„ + i ^fc„ + i ■ ( A7) 



In passing from Eq. (A4) to Eq. (A5), we have used the fact that the system is periodic in t, xi, .., a;„_i directions; 
the integral in Eq. (A5) is equal to Eq. (A4) because only the u-direction has boundaries. 

Since the 0's do not depend on the coordinates (a;„,-- - ,X2„), translation symmetry in these (n+1) directions 
remain intact, and the associated momenta p are good quantum numbers. As a result the 2n-dimensional system 

decouples into ni=n number of (n — l)-dimcnsional subsystems each labeled by a different p, and is described 

by the Hamiltonian H{k; m{p + 6{t, x))). Let S[n~i)D{p) be the topological action of each of these n — 1 dimensional 
systems (upon integrating out the fermions), they satisfy 



(A8) 



For a fixed 6{t, x) profile, the Pontryagin index of the map p — 5- h{'p + 9) = m{p + 9) / \rn{p + 9)\ is given by 
P= ' 



Area(5»+"T) J '^"^ V e^"-""+^«c, (p + 0)9^,71,, (p + 9) ■ ■ • ap„^,7i,„^,(p + 9), 
j^^^^ls^+i) J ^"^'P e'"'""^' K ip + 0)dg, n,, (p + 9)-.- 5,„^,n,„^, {p + 9)] 



(A9) 
(AlO) 



It is important to note that the Pontryagin index of rn{p + 9) is independent on 9. Now, we can rewrite Scs„ as 
follows: 



<S'cs„ = ScSr, 



? Area(Ui) / '^"^V — • ■ .5,„,,n,. 



(All) 



1 



Cn/P 



(n+ l)!(27r)" Area(S'"+i) 



■ 2n 



d"+'p / dtd^-'x / du 



^Ci---C„ + 2, 



s-^--'^^^e'''-'-^'daA,9a2eb2---da^+A^+, , (A12) 



E 



. Cn 



2tt 



P Area(S'"+i)(n+ 1) 



— j dtd"-^x J due''-'"+'e-'---^'n,,da,n,,---da^+,n,„+2 



(A13) 



In going from Eq. (A12) to Eq. (A13) we have applied the chain rule of differentiation. Note that because, for different 
p, the Hamiltonians H{k; Xrn{p + 9{t, x))) can be adiabatically deformed into one another without closing the energy 
gap, the subsystems correspond to different p should share the same form of WZW term with the same quantized 



17 



coefficient. Since 5cs„ = J2p^{n~i)Di'P); wc conclude that 



P Area(S'"+i) 



(A14) 



We note that in Eq. (A14) the dependence on p is impUcitly contained in h{p + 9{t,x)). By replacing m(p 
9{t,x,u)) with m{t,x,u), which smoothly interpolates between m(t,x,u = 0) = m(p) and rn(t,x,u = 1) = ■m{p 
0{t,x)), we obtain the following topological action upon integrating out the fermions in Eq. (1) with h{t,x,u) 
rfi(t, X, u)/|to(<, X, u)\ 



.Cn 27r '■^ 



S'(n-i)ri = -i^^^i^j^^^;;^ / rficf" / du e^^'^^^+^nc^dtric., ■ ■ ■ dun^^.^^^, (A15) 



which is exactly the WZW term we desire 
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